Mixed-state fidelity and quantum criticality at finite temperature 
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We extend to finite temperature the fidelity approach to quantum phase transitions (QPTs). This 
is done by resorting to the notion of mixed-state fidelity that allows one to compare two density 
matrices corresponding to two different thermal states. By exploiting the same concept we also 
propose a finite-temperature generalization of the Loschmidt echo. Explicit analytical expressions 
of these quantities are given for a class of quasi-free fermionic Hamiltonians. A numerical analysis 
is performed as well showing that the associated QPTs show their signatures in a finite range of 
temperatures. 
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I. INTRODUCTION 

Even, small and smooth changes in parameters gov- 
erning the dynamics of a physical system can, in some 
circumstances, result in a dramatic change of physical 
properties of the system itself. The traditional approach 
to these so-called critical phenomena is based on the no- 
tions of order parameter, correlation functions, symme- 
try breaking, and a general formulation in the frame- 
work of the Landau-Ginzuburg picture and renormaliz- 
tion group A phase transition can be triggered e.g., 
by a change of temperature of the system or, at zero tem- 
perature, by a change of some of the coupling constants 
( e.g., external fields) defining the system's Hamiltonian. 
In the first case one says that the transition is driven 
by thermal fluctuations, and the transition is referred to 
as classical, whereas in the second quantum fluctuations 
are held responsible for the transition and this latter is 
referred to as a quantum phase transition (QPT) (2j. 

In the last few years a big deal of interest has grown 
about the possibility of studying QPTs by means of ideas 
and tools borrowed from the new born field of Quantum 
Information Science [![. In this novel approach the key 
concept involved is quantum entanglement (or genuinely 
quantum correlations) and the idea is that quantum crit- 
icality can be suitably characterized in terms of the be- 
havior of different entanglement measures Q. Though 
the picture is still rapidly moving it is by now clear that, 
whereas there are no doubts that entanglement is indeed 
a valuable conceptual tool to analyze QPTs, one has, 
case by case, identify what is the entanglement measure 
e.g., block-entanglement or concurrence most suited to 
extract the relevant information. 

More recently an approach to QPTs based on another 
quantum information notion i.e., quantum fidelity, has 
been put forward [a, la] . The idea behind this novel ap- 
proach is quite simple: the dramatic chance of the struc- 
ture of the ground state occurring at the critical points 
can be fruitfully studied by analyzing their degree of dis- 
tinguishability. Since this latter quantity is related to the 



overlap i.e., scalar product, between two different ground 
states, the fidelity approach is basically nothing but than 
a metric one; the key ingredient being provided by the 
state-space distance between states corresponding to two 
slightly different sets of coupling constants. The expec- 
tation is that at the critical points a small change of the 
control parameters should result in an enhanced mod- 
ification of the state structure and this in turn should 
be detected by a greater state-space distance i.e., statis- 
tical distinguishabiliy, between the associated quantum 
states. In spite of its apparent naivety, this metric-based 
approach turns out to be able to provide an effective 
way to obtain, qualitative as well as quantitative, infor- 
mation about the zero-temperature phase diagram of a 
large class of non trivial quantum systems i.e., quasi- free 
fermionic systems @ and matrix-product states Q • The 
conceptual appealing of the fidelity approach for detect- 
ing boundaries between different phases lie in its univer- 
sal geometrical as well as information-theoretic nature; in 
principle no apriori knowledge of the symmetry-breaking 
mechanism and of the associated order parameters is re- 
quired. 

In this paper we are going to extend the fidelity ap- 
proach to finite temperature. This generalization can be 
achieved by exploiting the mixed-state fidelity introduced 
by Uhlmann Q and related to the statistical distance be- 
tween two density operators (Bures distance) . This finite 
temperature extension of the fidelity approach is a non- 
trivial technical step necessary to analyze the signatures 
of QPTs at non zero temperature and, more in general, to 
investigate the potential usefulness of the fidelity notion 
in the study of classical i.e., temperature-driven phase 
transitions 9]. In the following we will focus on the 
first task. Moreover, we will provide an explicit finite- 
temperature extension of another concept that has been 
recently used in the context of QPTs and ins pire d the 
whole fidelity programme: the Loschmidt echo [lOl [TTI| . 

The paper is organized as follows: In Sec. II, we an- 
alytically calculate the mixed-state fidelity of thermal 
states and a numerical analysis is used to demonstrate 



the signatures of QPTs at finite temperature. In Sec. 
Ill, we calculate the Loschmidt echo through a similar 
procedure. Sec. IV contains the conclusion. 

II. MIXED-STATE FIDELITY OF THERMAL 
STATES 

Let us consider the set of (mixed) quantum states 
S{H) := {p £ C(H) I p > 0,trp = 1}. The mixed-state 
fidelity is given by [8| 



F(po,pi)--=teypiPoPi- (i) 

This quantity measures the degree of distinguishability 
between the two quantum states p and p 1 . The fidelity 
is related to the statistical Bures distance: D(p ,p 1 ) = 
y/2(l — F). We will use Eq. (TT]) to compare two different 
thermal states 

P a = Z- 1 exp(-/3 Q ff Q ), 

Z a : = trexp(-(3 a H a ),(a = 0,l). (2) 

We define 

Fh , Hi (Po,Pi) :=F{Z^e-^o^ z -i e -^y (3) 
In particular one can consider the cases 

• (i) H a = H(X a ), \ 1 =\ + 6X; (3 = ft = 0. 
. (ii) H a = H(P a ),P 1 =/3 + S/3. 

The first case (i) is useful to study the finite- 
temperature signatures of a QPT occurring at some 
points, say A c , in the parameter space. The second 
case (ii) is considered for studying temperature driven 
phase transitions. Here the (3 dependence of the Hamil- 
tonian can be the result of the presence of e.g., a chemi- 
cal potential term (grand canonical ensemble) or of self- 
consistently determined coupling e.g., a mean- field BCS 
pairing term. In this paper we will focus on case (i), while 
the use of mixed-state fidelity for analyzing temperature 
driven PTs will be addressed in a forthcoming work. 

A. Commuting Hamiltonians 

Let us start by considering the simple case where the 
two Hamiltonians H$ and Hi commute; then 

Pi PoP? = PiPo = (ZoZ^- 1 eM-PoHo + (4) 
Therefore by using the definition {1} 

f (Pq,Pi) 
= (ZoZ^hrexp [-(p H + 
= (^ ^)- 1/2 ^exp(-(/3 OJ E;0+/3 1 ^)/2) ) (5) 




FIG. 1: (Color online) Mixed-state fidelity of finite- 
temperature thermal state of the XY model. The fidelity 
is a function of A and 7. Here, we choose the spin number 
N = 200 and the parameter perturbation S\ — Sj = 10 -2 . 
The temperatures of system in the above four figures are cho- 
sen to be j3 = 1, (5 = 10, f3 — 20, and (3 = 100, respectively. 
At low temperatures, e.g., (3 = 100, the mixed-state fidelity 
clearly shows the signature of QPT which occurs at absolute 
zero temperature. When the temperature increases, the de- 
cay of fidelity becomes less sharp and finally the signature of 
QPT disappears due to the thermal excitation. 



- V 












4 -ft* 













-1.5 -1 -0.5 0.5 1 1.5 



FIG. 2: (Color online) A cross section of Fig. 1 with 7 = 1, 
i.e., mixed state fidelity of finite-temperature thermal state 
of transverse Ising model. Here, the same as that in Fig. 1, 
the spin number is chosen to be N = 200 and the parameter 
perturbation 5\ = <?7 = 10~ 2 . Four curves represent four 
different temperatures: [3 — 1,(3 = 10, (3 — 20, and (3 — 100 
respectively. The decay of fidelity, though becomes less sharp 
with the increase of the temperature, well indicates the critical 
point of the transverse Ising model. 
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where H a \^f n ) = E" \^> n )(a = 0, 1). In particular if more- 
over H = Hi one immediately finds that fidelity of ther- 
mal states can be expressed entirely in terms of partition 
functions 



(0) 



We note in passing that this relation seems to suggest 
the possibility of making a direct connection between the 
metric/statistical notion of fidelity and purely thermody- 
namical quantities as well as the viability of the fidelity 
approach even for classical systems [9(. We would like 
also to observe that Eq ([6]) also gives the (pure-state) 
fidelity between the pure quantum states that one asso- 
ciates to the thermal states p and p 1 through the kind 
of classical-quantum correspondence discussed in [see 
Eq (1) there]. 

To further exemplify this commuting case let us con- 
sider diagonal fermionic Hamiltonians H a = ^2 k e^c^Cfe. 
One has 



pjfi) = z-\{ 

k 

= na+ 



From this one finds 



e -0a e k c h c k 



1 + (e 



-(/V°+/3i4)/2 



k J(l + e-W)(l+e-PW) 



(7) 



In the zero-temperature limit (3 — > oo, from the 
above equation, it is easy to check that Vk, e k e\ > 
=> Th ,Hi{oo,oo) = 1, wheras if 3k,e k e\ < 
Fh ,Hi (oo, oo) = 0. This situation is the one that encoun- 
ters in the AF-model analysis of [5] in the critical line 
with anisotropic parameter 7 = 0. There indeed when 
the magnetic field A (which plays the role of a chemi- 
cal potential in the fermionic picture) is changed in the 
range [—1,1] one has exactly that one of the single parti- 
cle eigenvalues changes sign and this results in a vanishing 
fidelity @. 



B. A more general case 

Now we move to consider a more general case directly 
relevant to the XY and mean-field BCS-like models. The 
Hamiltonian is given by 

H a = V H% = V e%(n k + n_ k ) + Aft-ic^ + h.c.), 



fc 



k 



where rife 



— J 



(8) 

J,Cfc and the c^'s are fermionic op- 
erators i.e., {cfc,c[,,} = 5k,k'- The Hilbert space 
factorizes H = <8>fcWfe. One has TL k ® 1~L-k = 
span{|00) k ,_ k ,|ll)k,-k,|01) k ,_ k ,|10) k ,_ k }. The first 



(last) two vectors span the even (odd) parity sector. The 
Hamiltonian ||5J) has a trivial action over the odd sector 
i.e., H k \ dd = e k lk,-k- Therefore, neglecting irrelevant 
constants one can write 



Hk — 2 e kJk + 2A£J fe , 



(9) 



where J? 



:= 1/2 (rife + n_ fe 
1/2 ^— ic£cL fe + h.c.J . For every k, the operators 

J k := l/2(e|,c_ fc + h.c), J k , J k span a su(2) Lie-algebra 
such that the even (odd) sector of 7ik ® 7~i-k is the 
J = 1/2 (J = 0) irreducible representation. Therefore, 
with obvious notation one can write 



hi = (e^ + A^)eo 2 = AKeo 2 



(10) 



where 



A 



\j ( e fe) 2 + (Afc) 2 , 



^k a k 



qa y 



9 a k : =tan- 1 (A^/e^). 

Moreover, 

expH^fr?) = ^(/3jeife,_fe, 

Pa = WiZtT 1 Q a k{P a )®H,-k (11) 
k 

where 

Qt{(3) : = cxp(-/3A£ J£) 

= cosh(/3A£) - J£ sinh(/?A£) (12) 

is a 2 x 2 operator in the even sector, and 

Z% = 2 + 2cosh(/3A£). (13) 

The fidelity for the two thermal states is then given by 

Fh ,Hi {fin, Pi) 

= II " Tfuff • ( 14 ) 

k V^k^k 

So, apart for the trivial terms in the odd sector in or- 
der to compute the fidelity, one has to consider the prod- 
ucts Q 1 k ((3- L /2)Q° k ((3 )Ql(f3 1 /2). As this is a product of 2x2 
matrices, one has 

= • v /Tr [^(/ 3 o)^(/?i)] + 2det(^(/3 )^.(/3 1 ) 

Note that here det(g k (f3 )gl(/3 1 ) = 1, then substituting 
the above equation to Eq. (fT4f leads to a simple expres- 
sion of the fidelity 



Fh ,hAPoi P\) 

2+ v /Tr[^(/3 )^(/3 1 )]+2 



n 



(15) 
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Now, we are left to compute the trace of Qk(fio)&k(fii)- 
The matrix product g k (P )g k (f3i) can be written as 



Tr [(PMaliPi)] 



Tr 



= Tr([cosh(/3 A°) - sinh(/3 A°K 

= 2 cosh(/3 A") cosh^A^) - sinh(/3 A" 

xTr ( e - iQfc<Tfc 



-PlAk a kz p-iakOk: 



(16) 



where 



The following trace can be evaluated as 



Tr 



-ia k CFki 



-Pl^kVkz p-iUk<Tk: 



= Tr ( e - i2a ^ [a kz cosh^A*) - sinn^A* )] ) 

= cosh(/3 1 A^)Tr(e" i2QfcCT ^cr fc2 ) - 2 sinh^A^) cos(2a fc ) 

= -2sinh(/? 1 A^,)cos(2a fe ) (17) 

Substituting Ea.([T7p to leads to 

= 2[cosh(/3 A«)cosh(/3 1 Ai) 

+ sinh(/3 A°) sinh^) cos(^ - 

Bringing all the terms together one eventually finds 

•^Ho.Hi (flo, Pi) 

= n{[ 1 + cosh (/9oA° fc )] [l + cosh^Ai.)]}^ 



x(l + -L[l + cosh(/3 A°) cosh^A 1 ) 
+ sinh(/3 1 A«)sinh(/3 1 Ai)cos(^-^)]^) 
It is easy to check that 



(18) 



^Ho,H! (13,(3) 



f}->c 



n 



' 1 + cos( 



n 



cos( - - fc ) 



(19) 



i.e., one recovers the zero-temperature result @. 



C. Numerical analysis 

We use XY model as an example to demonstrate our 
main idea. Here, two Hamiltonians are Hq = if (7, A) 
and Hi = ff (7 + £7, A + <5A) where 



H(i, A) = £ 



1+7 



L^a^ + Aaf 



(20) 



where 7 defines the anisotropy and A represents external 
magnetic field along the z axis. Obviously, when 7 = 1, 
the XY Hamiltonian ([2H)l reduce to the transverse Ising 
Hamiltonian. £7 and 5X represent small perturbation 
to the Hamiltonian. erf, a 6 {a;,y, z} are usual Pauli 
operators in the i-th lattice point. It is well known that 
this model, by means of a Jordan- Wigner transformation, 
can be mapped onto a quasi-free fermionic Hamiltonian 
of the type © [i.e., e k = cos - A, A k = 7 sin 2**]. 
Through a straightforward calculation, it can be obtained 
that A° k = A fe (7, A) and A\ = A fc (7 + £7, A + £A) where 



2kir 



+ 7 2 sin 2 



2kn 



(21) 



and 6% = cos- 1 [(cos ^ - A a )/A£], (a = 0,1, A = 
A, Ai = A + £A). We plot the mixed-state fidelity accord- 
ing to the analytical expression obtained above (fT8|) . We 
choose the number of spins to be 200, and the pertur- 
bation to be 10 -2 . In Fig. 1, we consider four different 
temperatures; = 1, = 10, = 20, and (3 — 100, re- 
spectively. At low temperature, e.g., (3 = 100, the sharp 
decay of fidelity in critical region clearly displays the 
QPT, which happens at zero temperature. What is more, 
the pattern of fidelity in low temperature is very similar 
to the ground-state fidelity [5j of the XY model. This 
similarity is natural because, as we mentioned above, 
the mixed-state fidelity approaches ground-state fidelity 
when the temperature decreases to zero. When the tem- 
perature increases higher, e.g., [3 = 10 or (3 = 20, the 
decay of fidelity, though less sharp due to the thermal 
excitation, still early shows the signature of QPT. How- 
ever, when the temperature becomes large enough, e.g., 
(3 = 1 , the dramatic decay of ground-state fidelity at crit- 
ical point is totally washed out by the thermal excitation 
and the signature of QPT at critical point eventually 
disappears. In Fig. 2, we plot the mixed-state fidelity 
of thermal state of the transverse Ising model. This is 
actually a cross section of Fig. 1 at 7 = 1. The decay 
of fidelity clearly indicate the critical point of transverse 
Ising system. 

It is important to stress that these numerical findings, 
besides illustrating the possibility of detecting and study- 
ing the finite-temperature signatures of QPTs, show that 
the ground-state fidelity approach to criticality is en- 
dowed with some robustness against "perturbations" of 
the ground state. In this case, in fact we have seen that 
mixing the ground state with excited eigenstates does not 
destroy the peculiar behavior of fidelity in the neighbor- 
hood of the QPTs i.e., the fidelity drop. 



III. LOSCHMIDT ECHO 

In Refs. [10( and [llj, the concept of Loschmidt echo 
has been used to investigate quantum criticality. The 
resulting picture is that the asymptotic value of the 
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FIG. 3: (Color online) Mixed-state Loschmidt echo finite- 
temperature thermal state of the XY model at an instant 
t — 10. Here, we choose the spin number N = 200 and the 
parameter perturbation S\ = Sj = 10~ 2 . The temperatures 
of system in the above four figures are chosen to be /3 — 1, (3 — 
10, j3 = 20, and f3 — 100, respectively. Similar to the mixed- 
state fidelity approach, the critical region is well indicated 
by Loschmidt echo at an instant when the system is at a 
low temperature. However, the signature of QPT disappears 
when the temperature increases to very high value. 



Loschmidt echo (for infinitely long times) is directly re- 
lated to the closeness of the system to the critical points 
i.e., the closer the systems to the QPT the smaller the 
asymptotic value of the Loschmidt echo. On the other 
hand, even the short time behavior brings about informa- 
tion about the QPTs. Indeed, for short times the decay of 
the Loschmidt echo appears to be gaussian ~ exp(— at 2 ), 
where the rate a has a diverging derivative (as a function 
of e.g., magnetic field) at the critical point fill ]. 

Now we show how to generalize the notion of 
Loschmidt echo in a natural way to the thermal state 
case. This natural extension to the realm of mixed-state 
can be obtained by using Eq. fTJ), 



exp(—iH a t) (a — 0, 1). 



In particular, one can chose p in — exp(— (3Hq) 
simplifies and one can define 



(22) 
then Eq. 



C Ho , Hl (P,t) = F 



Zo'F 



0H \Zo X ul{t)e- 0Ho U x {t) 
e-P Ha ,Ultt)e-P Ho Ui{tj 



— •^H ,ul{t)H U 1 {t)\^-'P)- ( 23 ) 

Then, from the above equation and Eq. I|15p . one finds 




FIG. 4: (Color online) A cross section of Fig. 3 with 7 = 
1, i.e., mixed-state Loschmidt echo (at an instant) of finite- 
temperature thermal state of transverse Ising model. Here, 
the same as that in Fig. 3, the spin number is chosen to be 
N = 200 and the parameter perturbation 5\ — Sj = 10 -2 . 
Four curves represent four different temperatures: = = 
10, P = 20, and /3 = 100, respectively. The critical point is 
well indicated by Loschmidt echo when the temperature is 
not too high. 



with 



e" 2 



<? k x p -t3A k o kz - 



Ux(t). (25) 



Notice now that the unitary operator Ui(t) can be writ- 



ten as U\(t) = e 



iAit 



Tr [qU(3)qI(P)] 

= Tr(e i ~^ <T ' cx e ~PA"a kz e -i^r-<r k a: e i^-<J k x e itAl<7 kz e ~-i -£-a kl 
°-<7kx p — /3A^cr fe;5 — i-f-cr kx p i-^-a kx — itA^a^ — i-§-cr hx 



xe 2 



Tr(e _ ^ A ^ CTfcz g iQ fc CT fex gi*Aj.<r fcz e -ia k a kl 
g-^A^o-fcj giQfcO-fcx e -itAl<j kz e ~ia k <J kx \ 

Tr(e^ A * CTfe e"^ A * vV ^ v ), 



(26) 



where V = e M "'"e"" A ''"'e"" ,t '' fc ". Using the formula 
e x = coshx — sinhir, one have 



Ch ,hAP, t)=U (ZlY 2 + \/Tr [QKPHifi)] + 2 



(24) 



Tr [qIWM] 

2 cosh 2 (j3A° k ) + sinh 2 (/?A°.)Tr(a fc2 V^a kz V) . (27) 



6 



The following trace is evaluated as 



Tv{a kz V^a kz V) 
Tr(cr fcz e iQfcCTfca =e 



itA. a k% „-ia k a k!c 



e 

Tr(cr fez e iQfe<T '" e itA * CTfc2! e ~ l2QfeCTfc * e ~ itK l°>~ e *«*<^ ^ ) 
Tr(e i2Qfc<Tfex e i * A k cr '" 



e ~i2a k a kx e ~itA k a kl 



Tr( [cos(iAi)+isin(A^)e l2Qfc,T ^(7 fcz e- l2afe ' Tfe -] 
x [cos(tAl)-isin(tAl)(7kz] ) 

2cos 2 (iAi)+sin 2 (<Ai)Tr(a fez e l2Q ^-a feze - l2Qfc —) 
2cos 2 (iA^)+2sin 2 (iA£)cos(4a fe ) (28) 



Substituting the above equation to Eq. (f27|) leads to 

= 2 cosh 2 (/3A°) + 2 sinh 2 (/?Ag ) 

x cos 2 (iAj,)+sin 2 (tAfc)cos(4a fe ) 

= 2 [1 - sin 2 (2a fe ) sin 2 (A^)] cosh(2/?A£) 

+2sin 2 (2a fe )sin 2 (A^). (29) 

Bringing all the terms together, we give an explicit 
expression for the Loschmidt echo for the class of Hamil- 
tonians HI 



£-Ho,fli (Pit) 

nti+ cosher 1 



(30) 



A' 

X{1 



V2 



{[l-sin 2 (0°-^)sm 2 (A^)] 



x cosh(2/3A°) + sin 2 (0° - B\) sin 2 (A^) + l} 1 / 2 } 

It is immediate to check that in the zero temperature 
limit one recovers the result obtained in [To| , 



C Ho , Hl (M H yj(l - sin 2 (#° - ^)sin 2 (A^). 



the increase of the temperature, the decay of the mixed- 
state Loschmidt echo at critical point becomes less sharp 
until finally disappears. Fig. 4 is a cross section of Fig. 
3 at 7 = 1, i.e., a mixed-state Loschmidt echo of trans- 
verse Ising model. Similarly to the mixed-state Fidelity 
of transverse Ising model in Fig. 2, the critical point 
is clearly indicated by the Loschmidt echo. With tem- 
perature increase, the decay of Loschmidt echo at criti- 
cal point becomes less evident, and finally disappears for 
sufficiently high temperature. 



IV. CONCLUSIONS 

In this paper we have shown how to extend to fi- 
nite temperature the fidelity approach to quantum phase 
transitions advocated in References [5[ and [6J] . This gen- 
eralization relies on the notion of mixed state fidelity ap- 
plied to (Gibbs) thermal states. Mixed-state fidelity is 
strictly related to the Bures metric measuring the statis- 
tical distance between two density operators, therefore 
this approach has a geometrical as well as an operational 
meaning and two of them are deeply intertwined. We 
provided an explicit analytical expression for both the 
fidelity and Loschmidt echo for an important class of 
quasi- free fermionic Hamiltonians including e.g., the XY 
model. A numerical analysis of these quantities has been 
performed; it clearly shows how the influence of the zero- 
temperature critical points extend over a finite range of 
temperatures. Besides representing a non-trivial general- 
ization of the former zero-temperature results, the find- 
ings reported in this paper suggest directions for further 
investigations. The most prominent being the possibility 
of using (mixed) state-space metrical quantities to study 
temperature driven phase transitions both in the quan- 
tum and in the classical case Q. The fact that simple 
geometrical notions seme to provide a unified frame for 
studying all these different kinds of critical phenomena 
is, we believe, conceptually quite appealing. 
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